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Abstract 

We show how an ion trap, configured for the coherent manipulation of external and internal 
quantum states, can be used to simulate the irreversible dynamics of a collective angular momentum 
model known as the Dicke model. In the special case of two ions, we show that entanglement is 
created in the coherently driven steady state with linear driving. For the case of more than two ions 
we calculate the entanglement between two ions in the steady state of the Dicke model by tracing 
over all the other ions. The entanglement in the steady state is a maximum for the parameter 
vàlues corresponding roughly to a bifurcation of a fixed point in the corresponding semiclassical 
dynamics. We conjecture that this is a general mechanism for entanglement creation in driven 
dissipative quantum systems. 

PACS numbers: 03.65.Yz, 42.50.Fx, 03.67.-a, 42.50.Vk, 
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INTRODUCTION 



Generically, many body quantum systems are known to be difficult to simulate efficiently 
on a classical computer. This is because the quantum system may explore regions of state 
space with non zero entanglement, giving these systems access to a vastly larger state space 
than is possible classically. In an open quantum system we may, in some circumstances, 
be able to resort to stochastic methods, such as Monte-Carlo simulations. However, this 
will not be possible for open systems in which the steady state itself is entangled, as in the 
example we describe here. Terhal et a/.|| have considered the possibility of using a quantum 
computer to simulate open quantum systems in thermal equilibrium. Plenio et al.^ have 
considered how decay can lead to entanglement rather than destroying it. Cabrillo et aZ.|| 
discuss creating entanglement in two or more àtoms, by driving the àtoms with a weak làser 
pulse and detecting the spontaneous emission. In a recent paper by Arnesen et al. |§], the 
authors look at a situation where the spins in a Heisenberg chain with an external magnètic 
field show entanglement in the thermal state with non-zero temperature. In this paper we 
consider the case of an ion trap quantum computer realisation to simulate the irreversible 
dynamics of iV two level systems. 

In the following we will investigate the use of an ion trap quantum computer as a simula- 
tion device for the so-called Dicke model J5|. This model includes resonance fluorescence of 
a set of two-level àtoms driven by a resonant coherent làser field as well as a collective decay 
mechanism. We will first describe how a collective decay mechanism may be realised for iV 
trapped ions interacting with a collective vibrational mode when the vibrational mode is 
subject to controlled heating. Simulating irreversible dynamics for a trapped ion has been 
previously suggested by a number of authors || . We consider the dynamics of the system 
relaxing from an initial state in which all ions are in the excited state. For two ions we ex- 
plicitly calculate the Wooters entanglement measure (concurrence) as a function of time, 
which shows that the system passes through entangled states on the way to the collective 
ground state. Next we add coherent driving to try to force the system into a non trivial 
steady state and show that, again for the case of two ions, this steady state can be partially 
entangled. Extending this result to many ions is not possible at present due to the lack 
of a general measure of entanglement of mixed states in higher dimensions. However, we 
can calculate the entanglement between two ions or àtoms by tracing over all the other ions 
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or àtoms. This will show us at least whether entanglement is present. Interestingly, the 
maximum entanglement occurs for parameter vàlues for which the corresponding semiclas- 
sical system undergoes a bifurcation and loss of stability of the fixed point. We conjecture 
that the loss of stability of a semiclassical fixed point will generically be associated with 
entanglement in the steady state of the full quantum system. 



THE MODEL 



In the seventies the Dicke model and cooperative effects were subjects of research in 
various groups (see e.g. [§, ^ [ÜJ and references therein). The model consists of a group 
of two-level àtoms which is placed in a volume with dimensions small compared to the 
wavelength associated with the atom's two-level dipole and evolves on time scales shorter 
than any J 2 breaking relaxation mechanism (see ||) like an angular momentum system 
which has collective atòmic raising and lowering operators J + and J_ with a fixed spin 
quantum number j = N/2, where N is the number of àtoms. 

In the rotating frame with Markov, elèctric dipole and rotating wave approximations and 
ignoring a small atòmic frequency shift the master equation for the density matrix of this 
group of àtoms under the cooperative influence of an electromagnètic field is 0, [ÏÏ], [Ï2[ 
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2J_pJ + — j + J_p — pj + J_ 



(1) 



where íl is the Rabi frequency and 7^ is the Einstein A-coefhcient of each atom. This model 



can be solved exactly [J IJ and it exhibits a critical-point nonequilibrium phase transition for 
Q/j = 7a in the limit íl, j — > 00 Q. 



Collective driving 

How do we get a similar master equation to Eq. ([!]) in an ion trap? The coherent evolution 
is easy: We just shine the same làser at the carrier frequency on all the ions at the same 
time, thus forcing each ion to undergo Rabi oscillations at the same frequency. If we start 
initially with all the ions in their electrònic ground state \g), the ions will not leave the 
j = N/2 space. From there we can then define collective angular momentum operators in 
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the following way 



jV 



(2) 



i=l 



N 



J. 



+ ' 



(3) 



í=i 



where the raising and lowering operators for each ion are defined by <r_ = |<7)(e| and à + = 
\e)(g\. With this the Hamiltonian for simultaneous resonant driving of all the ions can be 
written as 



where fi is the Rabi frequency for the electrònic transition. 
Cooperative damping 

For the collective decay mechanism we need to couple the ions equally to the same heat 
reservoir. In this paper we will argue that the reservoir may be taken to be the centre-of- 
mass vibrational mode. It is subject to heating and we assume that it is in a thermal state. 
To couple the ions to the vibrational mode we need another làser which, again, illuminates 
all the ions at the same time, but which is detuned from the carrier frequency to the red by 
the trap frequency so that the electrònic state of each atom gets coupled simultaneously to 
the centre-of-mass vibrational mode. This is described by a Hamiltonian for the ith ion of 
the form 



where we have introduced the bosonic annihilation operator a for the vibrational mode and 
the coupling constant is fi 2 = ^fio- The parameter rf = E r /(hMuj ) is the Lambe-Dicke 
parameter with E r the recoil kinetic energy of the atom, ujq is the trap vibrational frequency, 
and M is the effective mass for the centre-of-mass mode. The Lamb-Dicke limit assumes 
77 <§C 1 , which is easily achieved in practice. The frequency fi is the effective Rabi frequency 
for the electrònic transition involved. This sideband transition is used to efficiently remove 
thermal energy from the vibrational degree of freedom. If the rate of this cooling process can 
overcome heating due to external fluctuations in the trap potential, the ion may eventually 




(4) 




(5) 
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be prepared in the vibrational ground state. However, in general the vibrational state will 
reach a thermal mixture, p v = Z~ x exp[— fuvoà'à/kBT], where Z = Tr(exp[— ftwoà^à/kBT]) 
at some effective temperature T. If the heating and cooling rates are such that the system 
relaxes at a rate, a, large compared with any other time scale for ion motion, the ion can 
effectively be regarded as interacting with a thermal reservoir at temperature T. We can 
also arrange that the associated rate of energy dissipation is small, 7 a <C ujç, which simply 
requires that the coupling to the vibrational degree of freedom is weak. Finally we assume 
that the temperature of the vibrational degree of freedom is such that 7 a <C k B T/h. Under 
these assumptions we may eliminate the description of the vibrational motion from the 
dynamics and obtain a master equation for the electrònic state p^i of the ion 

^ = ~[H,p^\ +l A nV[^\P {Í) +1a(u + 1)X>[*IV° , (6) 

where the superoperator is defined by 

V[À] P = Àpiï - i(i f ip + ptiÀ) , (7) 

and where H is the Hamiltonian for any other reversible electrònic dynamics and n is the 
mean thermal occupation number of the vibrational degree of freedom. In what follows we 
assume that the cooling is very efficient and set n = 0. At any time we my turn off the 
cooling làsers, thus reducing 7a suddenly to zero. The master equation, Eq.(^), looks like 
the master equation for a single two level atom interacting with the many body radiation 
field via spontaneous emission and thermal absorption and we will discuss this case in the 
next section. In the case considered here, however, we ignore spontaneous emission which 
is justified if the electrònic states coupled by the làsers correspond to a long lived dipole 
forbidden transition. Thus the irreversible dynamics of the electrònic state is due entirely 
to the interaction with the phonons associated with the vibrational degree of freedom. 

If the external làser field on each ion is identical (in amplitude and phase) the interaction 
Hamiltonian is 

Éj = híl 2 (àJ + + a f J_) , (8) 

where J + and J_ are defined in Eqs. (|3]) and (0). For the case of a linear ion trap, with 
separately addressable ions, identical làser fields could easily be obtained by splitting the 
cooling làser into múltiple beams. In this way we can simulate an angular momentum system 



with quantum number j = N/2. This imposes a permutation symmetry on the system which 
reduces the effective Hilbert space dimension from 2^ to 2N + 1. Thus an exponentially 
large portion of the available Hilbert space, i.e. all the states with j < N/2, is not used 
in this simulation. However, it is easy to generate the relevant unitary transformations to 
simulate the j = N/2 angular momentum quantum system. 

It is not trivial to keep the vibrational mode in a thermal state of fixed temperature. One 
way of doing this was recently suggested by Kielpinski et al. [[y|: They put one ion, which 
is a different species of all the other ions, in the centre of a string of ions, so that they have 
an odd number of ions in the trap. Through this centre ion, which they can cool at will 
without disturbing the coherence of the other ions, the sympathetically cool all the ions and 
are thus able to keep the string of ions at a well defined temperature. The authors conclude 
that such a scheme of sympathetic cooling is "well within the reach of current experimental 



technique" [[nj. We assume for further calculations that the centre-of-mass mode is kept in 
such a thermal state by the outlined technique. In the following derivation of the master 
equation we do not explicitly put the làser cooling into the equation. We just assume that 
the vibrational state instantly, i.e. on a time scale fast compared to all the other processes 
involved, relaxes back into the thermal state. 

With all these assumptions we get the master equation describing the collective motion 
of the density matrix of all the ions 



dp _.n 

dt " ~~ ~ % 2 



J+ + J_,p 



+1A- (2J + pJ- - J-J + p - pJ-J^ 

•J~l — \ / „ * \ 

+ lA—(2J-pJ + -J+J-P-pJ + J-) , (9) 

where n is the mean phonon number of the vibrational centre-of-mass mode and 7a = 2ÇÍ%rf. 
Note that for n = this equation is identical to Eq. (|J) for the Dicke model. 



STEADY STATE AND ENTANGLEMENT FOR j = l 

With two ions we have j = 1 and from the master equation, Eq. (^), we can write down 
the equation of motion for the components of the three by three density matrix of the state 
of the system, taking into account that Tr(p) = 1 and that p is Hermitian. Getting the 
steady state is then a matter of simple àlgebra. 
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Once we have determined the steady state of the j = 1 system, we can rewrite this state 
in the underlying two-qubit basis. What we are interested in is the change of entanglement 
in the system as the parameters 7 and n change. The entanglement of two qubits is well 
defined [p], [TJ|, [T5| , and we choose the concurrence |7|] as a measure for it. 

A numèric evaluation of the concurrence leads to the plot in Fig. [TJ. What we see is 

cm 




FIG. 1: Plot of the concurrence as a measure of entanglement depending on the parameters 7 = 
7a/í^ and n, the mean phonon number of the thermal vibrational state. 

that we can get a certain amount of entanglement in the steady state of a coherently driven 
system that is coupled to a thermal reservoir. This is remarkable as the steady state is 
independent of the initial state, which can be unentangled. The coherent evolution alone 
does not lead to any entanglement for an initially unentangled state either, as it only consists 
of (simultaneously acting) single qubit rotations and no coupling between the qubits is 
present. Thus the entanglement is due to the cooperative decoherence in the system acting 
together with the coherent evolution. 

ZERO TEMPERATURE CASE. 

The analysis here is restricted to the case j = 1. For j > 1 and n^O, numerical methods 
will need to be employed to derive the steady state, however, in this case another problem 
will arise due to the fact that there is currently no measure of entanglement for N coupled 
qubits. Nevertheless, other phase transitions analogous to the one in the Dicke-model (see 
e.g. H |9| and references therein) will appear. 

For n = we can compare our results for the steady state to the ones calculated by Puri 
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and Lawande |L6| (see also Lawande et al. [fL7l). They calculate the steady state to be 



ps = è è (9*r m (9)- n J m Ji 

m,n=0 



(10) 



where 



2,k\9\ 



-2k 



(11) 



k=0 



is a normalization constant, g = i/y, where 7 = 771 /fi as defined above, and 

(2 + m + l)!(m!) 2 



H- 



2,m 



(2 -m)!(2m + 1)! ' 

With this we can write the density matrix of the steady state in matrix form as 



(12) 



Ps = D 



iV2y 



-2f 



\ 



iV2y 1 + 2 7 2 -iV2y - i2V2y 3 
\ -2 7 2 iV2j + i2V2y 3 1 + 2 7 2 + 4 7 4 j 

where we have calculated D as 

D = 3 + 4 7 2 + 4 7 4 . 



(13) 



(14) 



This gives the same result as above if we set n = 0. The concurrence for this special case is 

C{p) 0.12 r 




1 2 3 4 5 6 7 

1/7 1 7 



FIG. 2: Plot of the concurrence as a measure of entanglement depending on the parameters \g\ = 
I/7 = fi/ 7j 4 and l/\g\ = 7 (dashed line) with n = 0, i.e. for the Dicke-model. 

plotted in figure 0, but this time we plot it against \g\ = I/7 as well. Thus the Dicke-model 
as a special case of our model shows entanglement in the steady state. 
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For N > 2 up to date no definite measure of entanglement exists. But we can calculate 
the entanglement between just two of the N ions at least for the Dicke model, where the 
temperature of the bath is zero. The steady is then given by flTfl [LJl pf 



Ps 




^EE(-) 1^1 ■ ( 15 ) 

Z=0 Z'=0 

By writing this as a sum of states with angular momentum ji = 1 and 32 = j — 1 we can 
trace over the part of the Hilbert space with j 2 = j — 1 and thus get the density matrix in 
the steady state for just two ions (or àtoms in the original Dicke model). From there we 
can again calculate the concurrene. This time we plot it against the relative Rabi frequency 
|TÏ| Q r = We note that the màxima of the entanglement occur close to the critical 

point in the cooperative limit j, Q — > oo of the Dicke model, i.e. around íl r = 1. The two-ion 
entanglement is not the real measure of entanglement in the system. Thus we cannot take 
the cooperative limit as the two-ion entanglement goes to zero in this limit. However, we 
note from the plots in Fig. (||]) that the maximum value of the two-ion entanglement indeed 
does move closer to the point íl r = 1 for increasing N. 




0.5 1 1.5 

Q r = Q/(j 7 a) 



FIG. 3: Plot of the two-ion concurrence as a measure of entanglement depending on the parameters 
iïr = Í2/Ü7A) for (a) j = 1, (b) j = 4, (c) j = 16, and (d) j = 64. 



CONCLUSION 

In this paper we have demonstrated how the steady state of a dissipative many body 
system, driven far from equilibrium, may exhibit non zero quantum entanglement. This 
result is significant for two reasons. Firstly the steady state is a mixed state and the 



study of quantum entanglement for mixed states is a very active field of enquiry [18]. It 



immediately raises the question of whether the entanglement can be distilled and used 



as a resource for some quantum communication or computation task Secondly the 

maximum entanglement occurs at the same parameter vàlues for which the semiclassical 
dynamics of the system undergoes a bifurcation of the fixed point corresponding to the 
quantum steady state. At the bifurcation point the time constant associated with the fixed 
point goes to zero as the bifurcation is approached. This is reminiscent of a phenomenon 
that characterises quantum phase transitions, in which a morphological change in the ground 



state, as a parameter is varied, is associated with a frequency gap tending to zero p0|, pij . We 
conjecture that the association between the bifurcation of a fixed point of the semiclassical 
description and the maximum of entanglement will be a general feature of dissipative many 
body systems driven far from equilibrium. 
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